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ABSTRACT 


This report discusses the effects of the late-time high-altitude electromagnetic 
pulse (HEMP) on electrical transmission and distribution (T&D) systems. This 
environment, known as the magnetohydrodynamic electromagnetic pulse (MHD-EMP) is 
a very slowly varying electric field induced in the earth's surface, similar to the field 
induced by a, geomagnetic storm. It can result in the flow of a quasi-dc current in 
grounded power lines and in the subsequent magnetic saturation of transformers. This 
saturation, in turn, causes 60-Hz harmonic distortion and an increase in the reactive 
power required by generation facilities. This report analyzes and discusses these 
phenomena. The MHD-EMP environment is briefly discussed, and a simplified form of 
the earth-induced electric field is developed for use in a parametric study of transmission 
line responses. Various field coupling models are described, and calculated results for the 
responses of both transmission- and distribution-class power lines are presented. These 
calculated responses are compared with measurements of transformer operation under dc 
excitation to infer the MHD-EMP response of these power system components. It is 
found that the MHD- EMP environment would have a marked effect on a power system 
by inducing up to several hundreds of amperes of quasi-dc current on power lines. These 
currents will cause transformers to saturate which could result in excessive harmonic 
generation, voltage swings, and voltage suppression. The design of critical facilities which 
are required to operate during and after MHD-EMP events will have to be modified in 
order to mitigate the effects of these abnormal power system conditions. 
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MAGNETOHYDRODYNAMIC ELECTROMAGNETIC PULSE (MHD-EMP) 
INTERACTION WITH POWER TRANSMISSION 
AND DISTRIBUTION SYSTEMS 


1.0 INTRODUCTION 

A large-yield nuclear detonation at altitudes of several hundred kilometers above 
the earth will distort the earth's magnetic field and result in a time-varying geomagnetic 
field on the earth's surface. This varying geomagnetic field interacts with the finitely 
conducting earth to produce a time-varying electric field (E-field), also on the earth's 
surface. Known as the magnetohydrodynamic electromagnetic pulse (MHD-EMP), this 
E-field can induce a net voltage in long electrical conductors such as power transmission 
lines. If these conductors are electrically connected to the earth at both ends, a current 
can be induced in the conductors, causing damage or upset to certain electrical systems. 

This MHD-EMP environment, also referred to as the E 3 environment, is only one 
component of the total electromagnetic pulse, called the high-altitude electromagnetic 
pulse (HEMP), produced by a high-altitude nuclear detonation. The other components, 
denoted as E| and E 2 , have much larger amplitudes and a shorter duration. 
Consequently, they interact differently with the electrical power system. The MHD-EMP 
environment is characterized by a rather low E-field strength, on the order of several tens 
of V/km, and by a typical waveform pulse of several hundreds of seconds in duration. 
This part of the HEMP environment is similar to, but more intense than, that occurring 
for a solar geomagnetic storm. 

Very intense solar geomagnetic storms have previously upset long-line 
communications and power systems. It is predicted that the E 3 environment, which is 
much stronger than that of solar storms, will have significant effects. Most of the 
pievious solar storm effects have been noted for systems near the polar region, and some 
measures have been taken to mitigate these effects. The E 3 environment, however, could 
be experienced by power systems which are normally not affected by geomagnetic storms, 
and could therefore have an adverse effect on systems nationwide. Furthermore this E Q 
environment would be applied to the power system after the system has been excited by 
Ej and E 2 environments that have possibly been caused by other nuclear detonations. It 
is presently viewed that the synergistic effects of these environments will increase the 
likelihood of electrical power interruption in the event of a HEMP occurrence. 
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This report documents the results of a study on the interaction of MHD-EMP with 
power transmission and distribution (T&D) systems. Section 2.0 presents a brief 
overview of the MHD-EMP environments used in the study. Normalized MHD-EMP 
electric field waveforms are assumed to be simple excitation functions, which permits an 
analytical curve-fit representation for the time history of these fields. A more accurate 
description of this environment would involve a detailed transient waveform, its 
polarization, and its spatial dependence. 

t 

Section 3.0 describes the analytical models useful for coupling this MHD-EMP 
environment to T&D lines. Because of the quasi-static nature of the excitation, the 
models are essentially simple direct current (dc) circuit models. However, complications 
arise in attempting to model realistic line configurations having a large number of support 
towers and an overhead shield or neutral wires. These models are discussed in detail. 

Using the analysis models described in Section 3.0, a parametric study of a 500-kV 
power transmission line was conducted, and the results are documented in Section 4.0. In 
this study, the line length, the MHD-EMP waveform and amplitude, and the tower 
connections to the neutral conductor of the transmission line are varied. 

The possible effects of this HEMP environment on electric power systems are 
discussed in Section 5.0. Unlike the solar geomagnetic storm, which can last for many 
hours (or sometimes days), the E 3 environment for a single burst lasts for only several 
minutes. Consequently, direct damage to large power system components is not thought 
to be a problem. More serious, however, is the effect of power transformer saturation and 
the subsequent generation of 60-Hz harmonics. These phenomena can upset 
instrumentation in the system as well as lead to an increase in the reactive power required 
from the power generators. Both of these effects can lead to power system instability. 

Finally, Section 6.0 presents a summary of the findings of this report. The 
references are listed after Section 6. 
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2.0 THE MHD-EMP ENVIRONMENT 


2.1 Overview 


As discussed in ref.fl], the MHD-EMP or E 3 environment arises from a variation 
of the earth's magnetic field caused by a high-altitude nuclear detonation. The 
interaction of this time-varying magnetic flux field (B-field) with the imperfectly 
conducting ground causes a transient electric field to be induced on the surface of the 
earth, in a manner analogous to that occurring in a geomagnetic storm [2]. For a 
tangential B-field on the earth's surface, denoted by B(t), the corresponding electric field 
tangential to the earth is computed as 



(1) 


where a is the earth electrical conductivity and n = 4tt x io' 7 h/m is the permeability 
of the soil [1]. This resulting E-field is orthogonal to the B-field on the earth's surface 
and is much slower than the early-time E, HEMP environment. Typical waveform 
times for the MHD-EMP environment are on the order of several hundreds of seconds. 

The E 3 environment may be divided into two parts, based on the postulated 
mechanisms of production. The first part, for times between 0 < t < 10 sec, arises from 
the initial nuclear burst and its interaction with the earth's magnetic field. This part is 
referred to as the blast-wave component. At later times for 10 < t < 500 sec, a, second 
contribution to the geomagnetic field variation arises due to the late-time atmospheric 
heave. These production mechanisms are illustrated and discussed in ref.[l]. 

The environment used in refs.[l] and [3] for an EMP assessment of commercial 
power networks was based on early work reported in ref.[4]. This work involved 

performing a numerical simulation of the time development of the disturbed atmosphere 
using a code named MICE. 
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Recent refinements in the theory of MHD-EMP production by Austin Research 
Associates [5] have led to an alternate computer model for predicting the blast-wave 
environment. In addition, there are some preliminary results for the late-time heave 
component of the E 3 environment. This work has been used to develop an updated 
composite MHD-EMP environment, which can be used to estimate the behavior of 
induced currents in long transmission and distribution lines. 

According to ref.[5], a simplified way of viewing the early-time, blast-wave E 3 
environment is to consider a quasi-static problem in which a magnetic dipole moment at 
the burst point is used to represent a perturbation source for the geomagnetic field. This 
dipole is oriented in a direction opposing the earth's magnetic field. Below this dipole, at 
an altitude of about 110 km, is a conducting region, or patch, which is created by 
downward-streaming x-rays from the detonation. For our model, this x-ray patch is 
assumed to be perfectly conducting, with no penetration by the early-time magnetic field. 
However, the B-field from the dipole moment does reach under the patch by flowing out 
around the ends of the patch and then reconnecting in the region between the patch and 
the ground. In some of the MHD-EMP literature, this process is variously described as 
"propagation" or "diffraction," but since these are nominally high-frequency concepts and 
the MHD-EMP is quasi-static in nature, the use of these terms is avoided here. 

Under the x-ray patch, the blast-wave E-field is observed to be smaller than 
outside this shielded region. The field is oriented primarily in the west-east direction and 
does not appear to vary drastically with position. Outside the shielded region, the E-field 
appears to fall off with distance away from the burst, with the largest field occurring just 
outside the x-ray shield. 

For a burst near the magnetic north pole, these electric fields are also 
predominantly in the west-east direction. Although the field variation with position may 
be computed using the elementary dipole model representing the center of the magnetic 
bubble, our view of the E 3 environment assumes that the E-fields outside the x-ray shield 
are also in the west-east direction. 
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2.2 Simplified MHD-EMP Electric Fields 


For the purpose of illustrating the coupling of MHD-EMP environments to power 
systems, several different E-field waveforms may be used. Figures 1 and 2 present two 
different waveforms for the early-time blast component of the E-field on the earth's 
surface. The MHD-EMP E-field, denoted by E (t), is normalized by a factor E 
Figure 1 is typical for observation locations outside of the x-ray patch, while Figure 2 is 
for a location under the patch. The normalization factor, E max , depends on many 
factors, including the burst yield and other parameters, the exact observer location, and 
the earth's electrical conductivity. 

Figure 3 presents a typical late-time heave contribution to the E-field. As in the 
previous figures, this waveform is also normalized to unity. This component of the 
MHD-EMP environment is believed to be strongest directly under the burst, falling off 
rapidly as the observer moves away from ground zero. 

The total MHD-EMP E-field on the ground consists of some suitable combination 
of Figures 1-3, depending on the actual location of the observer. For some locations the 
late-time component of the environment will be very small, but for others it can be 
substantial. As an example of such a composite MHD-EMP waveform, Figure 4a presents 
a complete normalized waveform. Its frequency-domain spectral magnitude is shown in 
Figure 4b; note that most of the spectrum is located well below 1 Hz. For power systems, 
this implies that the MHD-EMP waveform appears as a quasi-dc signal, and that dc 
circuit modeling concepts will be appropriate for calculating system responses. 

For reasons mentioned previously, the MHD-EMP E-field vector direction depends 
on the burst location relative to the magnetic north pole. As in the case of solar 
geomagnetic storms, for burst locations near the north pole, much of the earth's surface 
experiencing the MHD-EMP environment sees a predominantly east-west E-field. 
Consequently, for this study, the direction of the E-field is assumed to be east-west. 
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Figure 1. Normalized Early-Time MHD-EMP E-Field Outside the X-Ray Patch. 



Figure 2. Normalized Early-Time MHD-EMP E-Field Under the X-Ray Patch. 
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Figure 4. Normalized Composite MHD-EMP E-Field. 


8 
















2.3 Curve-Fit Representation of the E-Field Environments 


For some calculational procedures, an explicit expression for the waveforms shown 
in Figures 1, 2 and 3 is useful. This expression may be obtained by fitting the waveform 
components to a suitable functional form using commercial curve-fitting software [6]. 

A suitable fit for the time-dependent E-field environment outside the x-ray shield, 
illustrated in Figure 1, is given by ref.[6] in the form of a power series expansion: 

E o( r )/ E max = a +br+cr 2 +dr 3 +er 4 -|-fr 5 +gr 6 -l-hr 7 -l-ir 8 -l-jr 9 -|-kr 10 (2) 

for 0 > t > 3.75 seconds. Here r is a shifted time variable given by r = (t - 0.4) 
seconds, and the following terms are computed in the curve-fitting process: 


Coefficient 


a 

b 

c 

d 

e 

f 


i 

j 

k 


Value 

-0.00066220 

-0.16634115 

-0.89240008 

0.547750685 

-2.28495005 

3.687696276 

-2.39828948 

0.788893721 

-0.13926466 

0.012495192 

-0.00044058 


Figure 5 illustrates this curve-fit (the solid line), along with the original waveform 
data for the normalized MHD-EMP waveform. Clearly, for shifted times r > 3.75, the 

curve-fit begins to deviate from the actual waveform, indicating that the fit should not be 
used in this region. 


9 







Figure 5. Plot of the Normalized Curve-Fit (Solid) and Waveform (Data Points) for the 
MHD-EMP E-Field Outside the X-Ray Patch. 


For the early-time component of the E 3 waveform for observation locations under 
the x-ray shield, the following curve-fit representation is possible: 

E 0 (r)/ E max = a+br+cr 2 +dr 3 +er 4 +fr 5 +gr 6 +hr 7 +ir 8 +j7- 9 +kr 10 (3) 

for 0 < r < 10 seconds, and again r is a shifted time variable given by r = (t - 0.4) 
seconds. For this expression, the following terms are defined: 
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Coefficient 


a 

b 

c 

d 

e 

f 


j 

k 


Value 

-0.15830209 

2.806631745 

-10.7123990 

11.35846963 

-5.53539948 

1.485628864 

-0.23541413 

0.022193314 

-0.00117987 

2.99420E-05 

-1.9980e-07 


Figure 6 illustrates this curve-fit (the solid line), along with the original waveform 
data for the early-time blast-wave component of the environment shown in Figure 2. 


FILE rill OAT SHirrtD BY -.4 SECONDS 
1 Eqn 216 y=(a+bx*cx2+dx3 +ex 4 + 4 x 5 +qx 6 + hx 7 +ix8*j x 3+k*10) r2=0.986562859 
a=-0.15830209 b=Z 8066317 c=-10.712399 0=11.35847 e=-5.5353995 1=1.4856289 
g=-0.23541413 h=0.022193311 i=-0.0011798676 j=2.9911583E-05 k=-1.9978131E-07 



Figure 6. Plot of the Normalized Curve-Fit (Solid) and Waveform (Data Points) for the 
Blast-Wave Component of the MHD-EMP E-Field Under the X-Ray Patch. 
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The late-time heave contribution to the E 3 E-field under the x-ray shield can be 
expressed as a ratio of polynomials as 


E„M/E 


max 


(a+br+cr^+dr^+er^+fr 5 ) 

(g t -t-hr+ir +jr+kr 1 ) 


(4) 


for 0 < t seconds, where the shifted time variable is now given by r = (t - 22) seconds. 
For this expression, the following terms are defined: 


Coefficient 


a 

b 

c 

d 

e 

f 


i 

j 

k 


Value 

-0.00541799 

-0.03056179 

0.031475266 

0.000660923 

-0.00051925 

-5.0690E-06 

2.61420E-06 

2.18120E-08 

-5.2680E-09 

-3.4780E-11 

3.73010E-12 


Figure 7 illustrates this curve-fit (the solid line), along with the original waveform 
data for the normalized late-time blast-wave component of the environment of Figure 3. 
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Figure 7. Plot of the Normalized Curve-Fit (Solid) and Waveform (Data Points) for the 
Late-Time Component of the MHD-EMP E-Field. 
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3.0 ANALYSIS MODELS FOR MHD-EMP RESPONSES 


The MHD-EMP electric field environment described in the previous section is 
capable of inducing currents in long sections of transmission or distribution lines in 
electrical power systems. Unlike the early-time E x HEMP environment, which can induce 
large currents in conductors that are not connected to the ground, the E 3 environment 
will induce currents in the lines only if the lines are connected to the earth at two or more 
points. This is due to the quasi-static nature of the MHD-EMP fields. In this section, 
models of several different transmission and distribution lines are examined and the 
MHD-EMP-induced currents computed. 

3.1 Isolated Line Section 

Single-phase electrical transmission or distribution lines are not usually connected 
to the earth through low impedance loads at each end of a long line, as this will short out 
the transmitted power. Consequently, the E 3 environment is not generally a problem for 
these lines. However, three-phase power systems with grounded Y transformers at each 
end can have this grounding configuration, as shown in Figure 8. Under normal 

operation, the 60-Hz power currents flowing along the line, Ip ^ and 1^, are all 120° out 
of phase so that the total return current through the earth or other neutral is 
approximately zero. The MHD-EMP earth-induced electric field, E Q , acts over the length 
L of the line to produce a net quasi-dc excitation voltage of 

V s = E 0 (t)L . (5) 

Although the peak value of the MHD-EMP E-field is relatively small compared 
with the 50-kV/m early-time Ej field, the line length L can be very long. Consequently, 
the excitation voltage can be on the order of several kV in certain cases. Of course, this 
value depends on the line length, the line orientation, and the environmental details. 
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a. Three-Phase Line and Transformer Configuration 
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b. Simple Equivalent Circuit 

Figure 8. Three-Phase Line and Equivalent Circuit for Computing MHD-EMP-Induced 
Current. 
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This MHD-EMP excitation of the line will induce a current to flow in the neutral 
conductors of the transformer, as indicated in the figure. This current, denoted by I , will 
divide and flow through each of the transformer windings, flow along each of the phase 
conductors, and return to the earth through the neutral of the second transformer. This 
current may be calculated by the equivalent circuit shown in Figure 8b, where 
represents the footing or grounding resistance of the neutral conductors of the transformer 
to the earth ground, R y is the parallel combination of the three transformer winding 
resistances, and R^ is the parallel combination of the three line conductor resistances. 
This latter quantity depends on the length of the section of line being considered. Using 
the term r^ to represent the total per-unit-length resistance (in ft/km) of the three 
phase conductors in parallel, this resistive component has the value R^ = r^L ft. 

Specific values for R^, R^, and R^ depend on the voltage class of the transmission 
or distribution line being considered, as well as on the location and construction of the 
transformer site. An efficient power system design attempts to minimize each of these 
resistances. Unfortunately, this minimization will lead to large MHD-EMP-induced 
currents. For example, overall circuit resistances can be on the order of several ohms to 
tens of ohms, with resulting currents of several hundreds of amps. 

A simple analysis of the circuit in Figure 8b results in the current response 

EL 

L = --2-- ■ (6) 

2(Rj+Ry) + fjL 

For short lines, the induced current is seen to be proportional to the line length. 
However, for long lines when the term r^L dominates, the current is seen to approach a 
constant value, independent of the line length. 

To obtain an indication of possible current responses for different line 
configurations, several line classes have been examined. The first case involves power 
transmission lines of the 138-, 345-, and 500-kV class. For these types of lines, an effort is 
made to keep the transformer neutral grounding resistance low, usually between 0.5 to 1.0 
ft. For this study, therefore, the grounding resistance for these lines is assumed to be 
Rf = 0.75 ft. Table 1 presents the other appropriate line and transformer resistances for 
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these line classes. Typical lengths of these transmission lines are indicated as the 
"Nominal L." 


Tabic 1 

Parameters for Transmission Line Analysis 


Line Class 

Nominal L 

Conductor 

r L 

R v 


(km) 

Data 

(fi/km) 

(«) 

500 kV,200 MVA 

100-500 

4x583,200 CM* 
18.7 ACSR* Bundle 

0.0083 

0.060 

345 kV,100 MVA 

100-200 

1,414,000 CM 
Expanded ACSR 

0.0135 

0.025 

138 kV,100 MVA 

50-150 

397,500 CM 
ACSR 

0.0540 

0.010 

* 

( CM = 

circular mil; 

ACSR = aluminum 

cable, steel 

reinforced) 


The second case considered is for subtransmission and distribution class lines 
operating at 12, 25, 34, and 69 kV, assuming that the transformer neutrals are connected 
to a substation ground mat at both ends of the line. This configuration provides a 
grounding resistance of « 1 0. Table 2 provides the length and resistance data for this 
case. The value for r L for the 12-kV line assumes that a 2/0 phase conductor wire is 
used. 


Table 2 

Parameters for Subtransmission and Distribution Line Analysis 


Line Class 

Nominal L 
(km) 

Conductor 

Data 

r L 

(D/km) 

R y 

(SI) 

69 

30-60 

4/0 Cu 

0.062 

0.04 

34 

15-30 

3/0 Cu 

0.072 

0.40 

25 

5-20 

3/0 Cu 

0.072 

0.30 

12 

2-15 

2/0 Cu 

0.091 

0.14 
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The third case considered is for subtransmission and distribution class lines 
operating at the same voltage levels as in Case 2, but with the transformer neutrals 
connected to a grounding rod on one end of the line and to the substation mat on the 
other end. The grounding rod typically provides a footing resistance of = 5.0 to 20.0 il 
For this study, this value has been taken as Rj « 10 fl The other data for this case are 
the same as in Table 2. 

For Case 1, involving the 138-, 345- and 500-kV transmission lines, Figure 9 
illustrates the behavior of the MHD-EMP-induced neutral current, shown as a function of 
the line length, L, for the three different line classes. These currents are normalized by 
the MHD-EMP E-field, E Q . Thus, the time-dependent induced current is given by the 
sample waveform of Figure 4a, multiplied by the proper factor L/E in Figure 9 for the 
line having a length specified by L. 

Figure 10 shows the induced neutral currents for the 12, 25, 34, and 69 kV line 
classes for Case 2 (i.e., the transformer neutrals grounded to the substation ground.) As 
in Case 1, the current limiting effect is apparent, especially in the case of the 12-kV line. 
Figure 11 illustrates the corresponding responses for the same lines, with the alternate 
grounding scheme (Case 3). Neutral currents in this case are significantly lower than in 
the previous case due to the higher grounding resistance provided by the grounding rod. 

Eq. (6) indicates that for very long lines the induced current is dependent only on 
the per-unit-length line resistance as I C /E Q = l/r L . This current limiting is apparent in 
some of the curves in Figures 9-11. For the transmission and distribution lines considered 
above, Table 3 summarized these peak MHD-EMP-induced currents. 

Table 3 

Maximum Normalized MHD-EMP-induced Current 


Voltage Class 

V E o 

(kV) 

(A-km/V) 

500 

120.5 

345 

74.1 

138 

18.5 

69 

16.1 

34 

13.8 

25 

13.8 

12 

11.0 
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Figure 9. MHD-EMP-Induced Current for Transmission Lines (Case 1). 










Figure 10. MHD-EMP-Induced Current for Subtransmission and Distribution Lines 
(Case 2 - Substation Grounding). 
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Figure 11. MHD-EMP-Induced Current for Subtransmission and Distribution Lines 
(Case 3 - Stake Grounding at One End and Substation Grounding at the Other). 


21 







3.2 Considerations for Shielded Lines 

A frequent practice with long transmission lines is to include one or more overhead 
shield conductors for lightning protection. Typically, these conductors are smaller than 
the power phase conductors; consequently, they have a larger per-unit-length resistance. 
These conductors are grounded at each tower, providing a conducting path to local 
ground. The MHD-EMP earth-induced E-field can also induce currents in these shield 
wires, and if the 3-phase power system is connected in some manner to the shield system, 
these currents can influence the level of current flowing in the transformer neutrals. 

Consider first the case of support towers and an overhead shield wire added to the 
3-phase line previously shown in Figure 8a. The overall line length L has N support 
towers, with a distance / between them. The phase conductors are supported by the 
towers, but do not have electrical connections to them. An overhead shield wire having a 
per-unit-length resistance of r g ft/km is connected electrically to the towers, and each 
tower is assumed to have a grounding resistance through the earth of R t ft. In this 
example, the transformer neutrals are assumed to be connected to the earth through a 
grounding resistance at a point other than the end support tower. This line 

configuration is illustrated in Figure 12a. 

The equivalent circuit for this line configuration is shown in Figure 12b. The N-l 
loops, or cells, in the part of the circuit representing the shielding line are caused by the 
sections of line between each of the towers. An induced voltage source of magnitude E Q / 
volts and a shield wire resistance of r / ft exist in each cell. The portion of the circuit 
representing the combined 3-phase conductors is the same as that of Figure 8b. Because 
the phase conductor line and the shield wire/tower line are assumed to be connected to 
the earth independently, there is no interaction between the two parts of the line. 
Consequently, the determination of the induced line current I proceeds as discussed 
earlier, using Eq. (6). 
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- L - 

a. Physical Line Configuration 


R Y R l = r L L = r L (N— 1 )l R y 



L-► 


b. Electrical Circuit Model 

Figure 12. Three-Phase Line with N Towers and Overhead Shield Wires Not 
Electrically Connected. 
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The decoupling of the phase conductor and shield line responses is based on the 
electrical isolation of the grounding resistances at the ends of the lines. In Figure 12a, the 
grounding resistances are shown as being completely isolated, a condition that is true only 
if the transformer ground is located far away from the end tower ground. When the 
transformer neutral ground is in close proximity to the end tower, however, there can be a 
mutual interaction between the two grounding resistances. This interaction may be 
modeled using a mutual resistance between the two parts of the circuit, as described by 
Sunde [7]. The effect on such mutual coupling will be to reduce the levels of the induced 
current in the transformer. Thus, considering the isolated (or unshielded) line provides a 
worst-case estimate of the response for these shielded lines. 

Another possible shielded line configuration is shown in Figure 13a, where the 
transformer neutral is connected to the tower grounds at locations A and A'. In this 
configuration, the 3-phase section of the line and the shield wire circuit are tightly 
coupled, and the transformer current response, I , will be greatly influenced by the towers 
and their grounding impedances. Figure 13b shows the equivalent circuit for this line. 
The connection between the transformer neutral and the tower ground is represented by 
the resistance element in this figure, although this resistance is generally very small 
and is usually neglected. 

In principle, the determination of the current I c is a straightforward task using 
circuit analysis. Unfortunately, for more than two of these towers, obtaining an analytic 
expression for the current is not feasible, due to the complexity of the circuit. As a result, 
the response must be obtained using numerical methods, and must be displayed 
parametrically. This may be accomplished by first isolating the effect of the source 
excitations, the towers, and the shield wire as an equivalent Thevenin circuit seen from 
the terminals A and A'. This is illustrated in Figure 14a. Once the open circuit voltage, 
V QC , and the input resistance, Rj n , are determined, the current I can be determined 
using Figure 14b as 


R in +2 ( R i +R y )+ r L L 


( 7 ) 
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Overhead Shield Wire 


Tower 1 


N — 1 


N 



a. Physical Line Configuration 


R Y R l = r L L = r L (N— 1 )l Ry 



b. Electrical Circuit Model 

Figure 13. Three-Phase Line with N Towers and Overhead Shield Wires with 
Interconnected Grounds. 
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a. Equivalent Circuit at Terminals A—A' 



b. Response Circuit for I 

V 

Figure 14. Analysis Procedure for the Current I c - 
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3.2.1 Solution for a Small Number of Towers 


The determination of and Rj n for the equivalent circuit of the tower shown in 

Figure 14a can be accomplished using mesh circuit analysis [8]. For N towers, this 

results in an NxN matrix that must be filled and inverted numerically. Two different 

forcing vectors are used, one for the determination of V and the other for Z. 

oc in 

For a moderate number of towers, say on the order of 100 or less, the solution of 
the circuit equations is possible with minimal computer resources such as a desk-top 
computer. Figure 15a shows the results of computing the normalized open circuit voltage 
V 0C /E 0 L 3S a function of the number of tower sections (i.e., the number of towers minus 
1). This voltage is shown for a range of different values of the dimensionless parametric 
ratio r g L/R t , which covers a wide variety of different circuit element values. As a result 
of the normalization, the maximum value of the open-circuit voltage response is unity, 
and the effects of the shield line and towers serve to reduce the effective voltage exciting 
the 3-phase conductors. Note that in the limiting case of the tower resistance R^ 
approaching zero, the parameter r g L/R^ becomes large, and the normalized excitation 
voltage is unity, with no effect of the overhead shield wires. 

The input impedance for the equivalent circuit in Figure 14a is shown in Figure 
15b, also in a normalized form as R in /R t - Input impedance is also plotted as a function 
of the number of tower sections, with the factor r g L/R t as a parameter. Note that for 
all cases the presence of this R jn impedance element in the response circuit of Figure 14b 
will tend to reduce the induced current response levels from that of the isolated 3-phase 
line. 


The presence of the overhead shield wire system is seen to affect the induced 
current flowing into the transformer neutrals in two ways: by decreasing the effective 
voltage source exciting this current, and by increasing the resistance seen at the source. 
In both cases, the overall effect is to reduce the net flow of MHD-EMP-induced current. 
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Figure 15. Behavior of the Equivalent Circuit at A and A' for Different Numbers of 
Tower Sections, Using Circuit Analysis. 
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The values of V Qc and Z^ n in Figure 15 can be used to compute the current 
induced in the transformer neutrals using Eq. (7). Figure 16 is a contour plot of the 
normalized current I^/E^ for the 138-kV line discussed in the previous section. The 
transformer and line resistance values used are those in Table 1, and an assumed tower 
footing resistance of R t = 25 fl was used. This plot shows possible current responses for 
any combination of line length and number of tower sections for line lengths 1000 km or 
less, or for 94 line sections or less. 

In this plot, several general trends are noted. For a line of fixed length, adding 
more grounded tower supports is seen to increase the coupled current response. This is 
due to the fact that as more tower impedance elements are added in parallel to the 
shielding part of the circuit, the effective input impedance, Z jn , is reduced. For the case 
of a fixed number of towers and an increasing line length, the effective excitation voltage 
increases due to the longer collection length of the line. Consequently, the induced 
current also increases. With any combination of length and number of towers, however, 
the induced current is lower than that which would exist without the shielding structure. 

3.2.2 Solution for a Large Number of Towers 

Practical transmission or distribution lines have support tower spacings that are 
relatively small compared with the overall line length. For example, the 138-kV line 
previously discussed has a typical tower span of 900 ft. (274.3 m); thus, for a 250-km line, 
approximately 910 towers would be present. Clearly the data presented in Figure 16 will 
not be adequate for predicting the line responses in this case. Moreover, if the circuit 
analysis approach were used to compute the V Qc and Z jn quantities directly, large 
computer resources would be required. To treat the problem of a 1000-km line with 
about 3650 towers, the solution of a 3650 * 3650 matrix equation would be required, 
which would pose practical problems. What is required is an alternate model for a large 
number of towers. 
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Figure 16. Contours of Normalized MHD-EMP-Induced Current I /E (A-km/V) for 

c o 

the Shielded 138-kV Line, Using the Circuit Model. 
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To conduct an analysis for a large number of towers, it is possible to consider the 
tower effects as being distributed along the length of the line, not as being individual 
lumped elements. This is done by defining a per-unit-length tower conductance of 


g t = NG t /L = N/(R t L) 


where N is the number of towers. This conductance may be used along with the 
per-unit-length resistance of the shield conductor r g to derive a set of coupled differential 
equations for the current and voltage distributions along the shield conductor. Figure 17 
shows this continuously loaded line running from terminal A at x = 0 to terminal A' at 
x — L. At an arbitrary point on the line, the line is considered to be made of the 
differential ladder circuit shown in this figure. By writing the Kirchoff voltage and 
current equations for this line section and taking the limit of dx -» 0, we derive the 
following set of differential equations: 


dX + r s I(x) = E q 


(9a) 


and 


dl + g t V(x) = 0 . 


(9b) 


To obtain a solution for V(x) at any point along the structure, we must first 
consider the solution of the homogeneous versions of these equations. These can be 
manipulated into uncoupled, second-order differential equations of the form 



(10a) 


and 



(10b) 


where the constant a has the value 



(11) 
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Figure 17. Continuous Model of the Overhead Shield Wire and Tower Supports. 
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The solutions to Eqs.(9a) and (9b) have the general form 


V(x) = ae- QX +be +ax 


(12a) 


and 


I(x) = — e -ox - — e~ ax 
R R 


(12b) 


where the constant R has dimensions of resistance and is defined as 




(13) 


The constants a and b are determined by the boundary conditions at x = 0 and L, 

and by the excitation term E . 

o 

As a practical matter, a realistic power line does not always have all towers spaced 
uniformly. However, this analysis is based on the averaged effects of the tower 
connections to the earth, and effects such as non-uniform spacing become unimportant as 
the number of towers becomes large. 

3.2.2.1 Determination of V 

oc 

A particularly useful solution to these equations is the Green's function [9], which 
is the response of the line due to a single excitation of the form E Q = V £(x-x ). This 
expression corresponds to a single voltage source of strength V g located at x = x . 
Once this generalized solution is obtained, the solution to an arbitrary distribution of 
excitation sources may be obtained by integration. Because the current on the shield at 
points A and A' is zero, it is possible to derive the following Green's function for the 
voltage response: 


e~ ax + e “( x - 2L > 



(14a) 


for the observation point x > x , and 

s 
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(14b) 


G(x;x g ) = • 2 V s 


<*(x-L) + e -a(x+L) 


(1 - e“ 2aL ) 


(e«( L - x s ) - e -a ( L-x s)) 


for x < x g . 

The excitation of the shield line running from 0 to L is due to the uniformly 
distributed voltage sources E Q dx. Using our Green's function, the solution for the 
complete line voltage at any point is expressed as 

L 

v (x) = J G(x;x g )E 0 dx g . (15) 

o 

This integral may be integrated analytically for G given in Eq. (14) to yield the 
following solution for the voltage: 


E„ 1 
V(x) = -5 


2a (1 — e -2oL ) 


[l + e -2a(L-x)j r (1< -ax } + e -2ax ,^+ax. 


- (l + e —2ax ] + e- 2 ^^) (l-e +a < L - x >)] 


Evaluating this expression at the end points of the line at A' and A gives 


(16) 


and 


V(L) = 


E ° (1 ~ e ~ aL ) 2 
a(l - e -2orL ) 


V(0) = - V(L) . 


(17a) 

(17b) 


These are the line-to-earth voltages. What is needed for the Thevenin circuit 

voltage is the voltage between these points, which is the V quantity shown in Figure 

oc 

14a. By applying Kirchoff's voltage law to the loop A-O-L-A' in Figure 17, the open 
circuit voltage may be expressed as 
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or 


V oc = E o L + v (°) - V(L) 


(18a) 



= E 0 L 


2 (1 - e -aL ) 2 
aL (1 - e -2oL ) 


(18b) 


3.2.2.2 Determination of R. 

in 

The input resistance between terminals A and A' is also needed for this analysis. 
At these terminals, the line of Figure 17 can be represented by the generalized 2-port 
circuit shown in Figure 18. 


A 

1(0)-t* 

V(0) 

4 -, 

i Zj 

A’ 

-t'(l) 

V(L) 

A 

0 


L 


Figure 18. Two-Port Representation of the Shield Line Between Points A and A'. 


This circuit can be represented by the open circuit impedance parameter matrix 
[10] of the form 


r-H 

> 

— 

Z 11 z 12 


V 

V 2 


Z 21 z 22 


>2 


where z 12 z 2 j by reciprocity, and z^ = z 22 by structural symmetry. The 

impedance elements in Figure 18 are related to the impedance parameters in Eq. (19) by 
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and 


Z 1 - Z 2 “ Z 11 ’ z 12 

(20a) 

Z 3 = z 12 . 

(20b) 


The net resistance between A and A' when both of their terminals are open circuited is 
immediately seen to be 



“ Z 1 + z 2 - 2 ( z n * z i 2 ) • 


( 21 ) 


The impedance parameters and z^ for the line may be calculated by placing a 
voltage source at terminals A-0 of Figure 18 and computing the input current Ij and 
the transfer voltage V 2 with the open circuit condition ^ = 0 at terminal A'-L. Using 
Eq. (12) with the end condition that I(L) = 0, the following solutions for V(x) and I(x) 
are obtained for. the case of a lumped voltage excitation of Vj at x = 0: 

V(»)- Vl -i ta u («~” + «~ aaL «”) (22a) 

(1 - e zaL ) 
and 

I(x)--- • (22b) 

fl(l - e zaL ) 

Thus, and z^ may be evaluated as follows; 


and 


V(0) 

1 ( 0 ) 


I(L)=0 


(1 + e~ 2aL ) 
% -e- 2aL ) 


z 12 ~ Z 21 ~ 


V(L) 

1 ( 0 ) 


I(L)=0 


2 R e aL 
- e _2aL ) ' 


(23a) 


(23b) 


Combining Eqs. (23a) and (23b) with Eq.(21) thus gives the following expression for the 
input resistance: 



e-“ L ) 2 

r 5 ®)' 


(24) 
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3.2.3 Sample Line Responses 


Equations (18b) and (24) are the Thevenin equivalent circuit parameters for the 
circuit between A and A' for the case of a large number of towers. For small values of N 
these results should agree with the results of the lumped circuit analysis discussed in the 
previous section. Both of the parameters aL and R in Eqs. (18b) and (24) may be 
written in terms of the dimensionless parameter r g L/R t used in Figure 15 as 


and 


ah = 


N 


V; 

R. 


(25) 


Figure 19a presents the normalized quantity V oc /E q L from Eq. (18b) as a function of 
the number of line sections (N-l) with the factor r g L/R t as a parameter. These data 
compare very well with the discrete-line model results in Figure 15a. Figure 19b presents 
the input resistance data for the same case using Eq. (24). For low values of the factor 
r s L / R t , these plots agree very well with those in Figure 15b. However, as this factor 
increases, there is a slight difference between the discrete and the continuous models. 


As another check of the validity of the formulation for the continuous line model, 
the calculation of the current in the 138-kV line discussed previously is presented in 
Figure 20 for up to 94 tower sections. This calculation compares very well with that in 
Figure 16, showing the same family of responses using the discrete circuit model. Note 
the differences in the curves occurring for small values (i.e., on the order of 10 to 20) of 
tower sections where the continuous model is not valid. 
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Figure 19. Behavior of the Equivalent Circuit at A and A' for Different Numbers of 
Tower Sections, Using the Continuous Model. 
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By way of comparing calculation times for the results in Figures 16 and 20, the 
lumped circuit model required about 12 hours of CPU time on an 80386/80387 computer 
running at 33 MHz. The distributed circuit model, however, required only about 30 
seconds for the same set of calculations, and this time does not increase as the number of 
tower sections increases. Therefore, this model is useful for treating the more realistic 
case of several thousand towers on the line. 

To provide a further indication of the MHD-EMP responses for shielded lines, a set 
of calculations for the 7 different voltage classes of transmission and distribution lines 
analyzed previously have been considered. Table 4 presents typical tower span distances 
for these lines, the assumed tower footing resistance (R t ), the longest length of line and 
the number of towers in that length (longest L and N), and the resulting per-unit-length 
tower conductance parameter (g t ). 


Table 4 


Parameters for Shielded Line Analysis 


Voltage Class 

Span 

R t 

Longest L 

N 

h 

(kV) 

(m) 

(ft) 

(km) 


(mhos/km) 

12 

60.90 

50 

100 

1642 

0.328 

25 

76.20 

50 

100 

1312 

0.262 

34 

91.40 

50 

100 

1094 

0.219 

69 

152.4 

50 

100 

656 

0.131 

138 

274.3 

25 

1000 

3645 

0.146 

345 

304.8 

25 

1000 

3280 

0.131 

500 

304.8 

25 

1000 

3280 

0.131 

Using the 

distributed tower model. 

, we have calculated the MHD-EMP-induced 

current in the above 7 classes 

of shielded lines. Tables 

1, 2, and 

4 summarize the 

parameters used 

for each line. 

Figures 

21a-g show the 

contours 

of the normalized 


transformer neutral current I C /E Q for various combinations of line length and number of 
tower sections. The straight line in each of these figures represents the locus of points 
corresponding to the tower span lengths given in Table 4 above. Other span lengths 
would be represented by similar straight lines having different slopes. 
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a. 500-kV Class Line 



b. 345-kV Class Line 

Figuie 21. Contours of Normalized MHD-EMP-Induced Current I /E (A*km/V) for 
Lines with Grounded Shield Wires. 
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Figure 21. Contours of Normalized MHD-EMP-Induced Current I C /E Q (A-km/V) for 
Lines with Grounded Shield Wires (Con't). 
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f. 25-kV Class Line 


Figure 21. Contours of Normalized MHD-EMP-Induced Current I /E (A-km/V) for 
Lines with Periodically Grounded Neutrals (Con't). 
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Figure 21. Contours of Normalized MHD-EMP-Induced Current I /E (A-km/V) for 

V U 

Lines with Periodically Grounded Neutrals (Concluded). 
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3.3 Induced Current on Grounded Neutral Conductors 


Judging from the dc circuit models developed in the previous sections, we can see 
one way of mitigating the effects of the MHD-EMP environment on the power system is 
to use an ungrounded A transformer at one end of the transmission line, as shown in 
Figure 22a. Because there is no dc path to ground through this transformer, there will be 
no dc current flow through the transformer windings. Consequently, no adverse effects on 
the transformer will be noted. 


An MHD-EMP electric field, however, can interact with the periodically grounded 
neutral conductor to produce a net current in the neutral conductor. This current is 
denoted by I n in Figure 22a. If the neutral conductor is connected to the ground far 
from the facility at the end of the line, there will be no adverse effect on the operation of 
internal power equipment. However, if the neutral conductor penetrates into the facility, 
there will be an injection of the MHD-EMP current into the facility. The effect of this 
current on the equipment within the facility will depend on the electrical connectivity of 
the ground conductor and the other internal equipment. Consequently, the effect cannot 
be predicted without detailed system information. 

It is useful to predict the possible levels of the MHD-EMP current which could be 
injected into a facility by the grounded neutral conductor. Figure 22b shows an 
equivalent circuit useful for describing the current. This circuit is similar to that in 
Figure 12b except it has the phase conductors removed and two additional resistances 
added at each end of the periodically grounded line. These resistances are denoted by R f 
and represent the footing resistance of the neutral conductors to the earth. The other 
circuit parameters are the same as in Figure 12b. 

The determination of the current flowing in the neutral line to ground, I , may be 
evaluated either by using dc circuit analysis or by using the analytical model discussed in 
the previous section. As in the previous examples, the circuit analysis is feasible only for a 
small number of tower sections, with the analytical solution being better suited for the 
case of many towers. 
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a. Physical Line Configuration 




b. Electrical Circuit Model 

Figure 22. Three-Phase Line with Un-Grounded A Transformer and Periodically 
Grounded Neutral Conductor. 
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Using the expressions for the voltage and current on the overhead neutral wire 
given in Eqs. (12a) and (12b), along with the loading of the line by Rj at each end, the 
following analytical expression for the neutral current, I , can be developed: 

E 0 L(l-e- aL ) 2 

0 r s L(l-e _aL ) 2 + aLR t (l-e -2aL ) (27) 

where the parameter aL is given by Eq. (25). 

There are 6 independent parameters for evaluating this expression: E q , L, R f , R , 
r s’ and N. To generalize this result slightly, we can plot the current in the normalized 
form I n /E 0 as contours in the L-N plane, as was done in Figures 20 and 21. Doing this 
requires that the remaining 3 parameters be defined. The tower footing resistance is 
taken to be the same as in the previous study, i.e. R t = 25 fi. Similarly, the grounding 
resistance at each end of the line is assumed to be that of a good substation ground, i.e. 

Rf = 0 75 The neutral conductor is assumed to be a #2 stranded aluminum wire with 
a per-unit-length resistance (r g ) of 0.876 D/km. 

Figure 23 presents the behavior of the normalized neutral current, as calculated 
using Eq. (27). As a check of the calculation, the same case was run using the dc circuit 
analysis, giving virtually identical results. Assuming that this neutral line is similar to 
that of a 12-kV line having a nominal tower separation distance of 60 m, the straight line 
on the figure represents the locus of points having this tower separation. 

The data in Figure 23 can be viewed in several different ways to give insight into 
the behavior of the induced current. By taking vertical cuts along the lines of constant 
line length, a family of curves of the normalized current vs. the number of tower sections 
may be generated. These are shown in Figure 24. 

Figure 25 presents the normalized MHD-EMP current vs. the line length for a 
different number of tower sections. It is clear that the MHD-EMP current at the end of 
the line is reduced when the number of tower sections is increased, because of the 
shunting effect of having a large number of towers. 


47 







Figure 23. Contours of Normalized MHD-EMP-Induced Current I n /E (A-km/V) at the 
End of the Periodically Grounded Neutral Conductor. 
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Figure 24. Normalized MHD-EMP-Induced Current vs. Number of Tower Sections for 
Different Line Lengths. 




Figure 25. Normalized MHD-EMP-Induced Current vs. Line Length for Different 
Numbers of Tower Sections. 
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It is not realistic to increase the tower density on the line too much, however, since 
there is usually a minimum tower spacing distance permissible for each line class. Figure 
26 shows the normalized MHD-EMP current vs. the line length for several different 
section spacings. The curve marked "60" represents the result for the nominal tower 
spacing distance of 60 m, while that marked "L" represents the case of only two support 
towers at each end of the line. There is a reduction of about 30% in the MHD-EMP 
current induced in the neutral conductor in this case. 



Figure 26. Normalized MHD-EMP-Induced Current vs. Line Length for Different Section 
Spacings. 


3.4 Cumulative Probability Distribution for Line Responses 

The induced current responses discussed above are all based on the assumption 
that the line is oriented in the direction of the maximum E-field, that is to say, in the 
west-east direction. In any real case, the line could have any direction, and consequently, 
the induced transformer neutral current would be reduced by a factor cos(</>) where 4> is 
the angle between the field and the line. 

A useful way of describing this effect is to use a cumulative probability 
distribution [11] which presents the probability that a particular response will exceed a 
specified value. Figure 27 illustrates this function for the induced peak current, 
normalized to a maximum value of unity, for the cos(^) field variation. 
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This figure may be used to provide a rough indication of the probabilities of 
occurrence of an MHD-EMP response for lines either under or outside the x-ray patch. 
As previously mentioned, this model neglects the changes in strength of the MHD-EMP 
E-field outside the shield, as well as the changes in direction of the field. A similar 
cumulative probability distribution curve could be developed using data from a more 
detailed MHD-EMP environment calculational model if the E-field environments could be 
provided at a sufficiently large number of observation locations on the earth. 
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Figure 27. Universal Cumulative Probability Curve for Induced MHD-EMP Current 
Responses. 
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4.0 PARAMETRIC STUDY OF MHD-EMP EFFECTS ON A POWER 

TRANSMISSION SYSTEM 

4.1 Introduction 

To better illustrate the effects of the MHD-EMP environment on electrical power 
systems, a parametric analysis [12] was performed on a specific 500-kV power 
transmission line, using the estimated E 3 environments and the coupling models described 
in the previous sections. The behavior of the unloaded transformer magnetizing current 
was analyzed; the results are presented parametrically with respect to the MHD-EMP 
waveform type and amplitude and the transmission line total length and span length. 
The system configuration, the analysis overview, and the results of the parametric 
analysis are all described in more detail in this section. 

4.2 Description of the 500-kV Transmission System 

The transmission system was based on the Minnesota Power Company 500-kV line 
between Minneapolis and Dorsey, Minnesota. Specifically, the system consists of a 
500-kV transmission line, terminated by three-phase transformer banks at both ends. 
Figure 28a illustrates a single line diagram of the test system. Pertinent data for this 
transmission line are listed in Table 5. During normal operation, the sky wires are not 
electrically connected to the transmission line towers. Tower configuration data specifies 
the location of the center of each phase bundle and each sky wire with respect to a 
Cartesian coordinate system with its origin located at the center of the tower base. 

Each of the three-phase transformer banks consists of three single-phase 
transformers connected A/grounded Y. The grounded Y side is connected on the 500-kV 
transmission line (see Figure 28). The characteristics of each single-phase transformer are 
as follows: 


Voltage Rating 
Power Rating 


Leakage Reactance 
Magnetizing Current 
Winding Resistance 


115/288 kV 
350 MV A 
0.10 pu 
0.002 pu 
i.5 n 
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Table 5 

Data for 500 kV Transmission Line 


General Data 


Configuration 
Bundle Spacing 
Earth Resistivity 


3 conductor bundles per phase 
18 inches 
100 fl -m 


Conductor Data 


Conductor 

Type 

O.D. 

(inches) 

Resistance 

(f)/km) 

Ground 

Conductors 

7/16 steel 

0.4375 

2.76 

Phase 

Conductors 

1192 ACSR 1.3020 

Tower Configuration Data 

0.05 

Conductor 

x-Coordinate 

(feet) 

y-Coordinate 

(feet) 

Phase A 

-32.0 


97.5 

Phase B 

0.0 


97.5 

Phase C 

32.0 


97.5 

Sky Wire 1 

-35.0 


129.5 

Sky Wire 2 

35.0 


129.5 


(Sky wire = transmission-line shield wire) 
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500-kV line 




473 miles long 



A 


a. Single-Line Diagram 


500 kV line a 



b. Three-Phase Diagram 

Figure 28. Configuration of 500-kV Transmission-Line for Parametric Study. 


The transformer core magnetization characteristics are described by specifying the 
flux linkage as a function of the magnetization current. This function is illustrated in 
Figure 29. The numerical values of this curve have been provided by GE. 

The electric load of the system was neglected in the parametric analysis. 
Experimental evidence obtained from transformer testing in Minnesota [15] indicates a 
substantial effect on transformer saturation constants due to the electric load. This 
effect, however, is not included in this parametric study, but is discussed in Section 5.0. 
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Figure 29. Transformer Core Magnetization Characteristics. 
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4.3 Description of the Parametric Analysis 


The behavior of the 500-kV transmission system was analyzed in the presence of 
the E 3 excitation. A parametric analysis was performed with respect to the excitation 
E-field waveform type and amplitude, the transmission line total length, and the span 
length. 


The two E 3 waveforms, defined as early time and late time waveforms and shown 
in Figures 2 and 3, were used in this analysis. Both of these waveforms are shown 
normalized, so that their peak values are 1 V/km. The transmission line system was 
analyzed for several different assumed peak amplitude values for each type of waveform. 
Specifically, the waveform amplitudes shown in Table 6 were used. 


Table 6 

Peak MHD-EMP E-Fields Used in Parametric Study 


Early-Time Waveform 
Amplitude (V/km) 

10 

20 

60 

120 

300 

500 


Late-Time Waveform 
Amplitude (V/km) 

10 

20 

40 

80 

180 

300 


For each of these above E 3 waveforms, a parametric analysis was performed with 
the values of earth resistivity, tower footing resistance, transmission line total length, and 
span length shown in Table 7. 
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Table 7 


Line and Ground Constants Used in Parametric Study 


Soil Resistivity 
Tower Footing Resistance 
Line Length 
Span Length 


30, 50, 100, 200, 400, and 1000 km 
0.16 and 0.40 km 


100 fi-m 
25 O and » 


The effects of MHD-EMP on the test system are assessed by computing the system 
transient response. The results of this study are described in the next section. 


4.4 Parametric Study Results 

The power system transient response was computed for each combination of the 
selected parameters given in Tables 6 and 7. The resulting transformer magnetization 
current and flux linkage, as well as the transmission-line zero-sequence current, were 

computed as a function of time for each case. Examples of these results are plotted in 
Figures 30 and 31. 

Since the parametric analysis results consist of a large number of plots, only 
selected characteristics of the obtained plots were tabulated so that the results could be 
presented in a compact form. Specifically, the maximum transformer dc offset of the 
magnetization current (for each set of selected parameters) and the time delay to reach 
magnetic core saturation were tabulated as a function of MHD-EMP waveform type and 
amplitude, soil resistivity, transmission line total length, and span length. The dc offset 
of the magnetization current is a very important quantity, since it determines the 
magnitude of the harmonics generated by the transformer. 
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a. Transmission-Line Zero-Sequence Current 


o 



b. Transformer Magnetization Current 


o 



c. Transformer Magnetic Flux 


Figure 30. Transmission Line Response for the Early-Time E 3 Waveform (E 3 Amplitude 

= 10 V/km, p = 100 fl*m, Shield Is Grounded, Line Length = 30 km, and Span Length = 
0.16 km). 
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a. 1 ransmission-Line Zero-Sequence Current 


«o 



b. Transformer Magnetization Current 



c. Transformer Magnetic Flux 

Figure 31. Test System Response for the Late Time E 3 Waveform (E 3 Amplitude = 10 
= 100 ^ m ’ Shidd IS Grounded ’ Line Len g fch = !000 km, and Span Length = 
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Figure 32 illustrates the maximum dc offset of the magnetization current in 
amperes for MHD-EMP excitation using the early-time waveform type as a function of 
the transmission line length and E-field amplitude. Parts "a" and "b" of the figure show 
the results obtained assuming that the shield wires are grounded at all towers, which were 
located with span lengths of 0.16 km and 0.40 km, for parts "a" and "b," respectively. 
Figure 32 shows the results obtained assuming the shield wires are insulated from the 
supporting towers. 

For most parameter combinations and with early-time waveform excitation, the 
transformer did not reach saturation. Therefore, the time delay to saturation was not 
tabulated for the early-time waveform cases. 



a. Shield Wire Connected to Tower with Span of 0.16 km 

Figure 32. The Maximum Magnetization Current for the Early-Time MHD-EMP 
Waveform when p = 100 fi-m. 
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c. Shield Wire Insulated from Tower 

Waveform when*p “woTm. Ma « netization Curre “‘ fo ‘ ‘he Early-Time MHD-EMP 
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Figures 33 illustrates the maximum magnetization current (in amperes) flowing 
into the transformer for the late-time MHD-EMP waveform for various transmission line 
lengths and waveform amplitudes. Parts "a" and "b" of the figure present results 
obtained for the two tower spacings with the shield wire connected to the tower, and part 
"c" shows the results for the shield wires insulated from the supporting towers. 

The time delay from the initiation of the late-time MHD-EMP waveform to 
saturation was also calculated and plotted for each combination of selected parameters. 
Specifically, the time delay to saturation is defined as the time interval from the 
initiation of the excitation waveform to the time instant that the transformer 
magnetization current reaches one tenth of its maximum value. The time delay to 
saturation as a function of transmission line length and the assumed E-field waveform 
amplitude are shown in Figure 34. The indicated time delay is in seconds. For some 
cases involving the smaller MHD-EMP amplitudes, transformer saturation was not 
obtained and the plots are truncated. Parts "a" and "b" of Figure 34 are for the shield 
wires grounded at every tower, and part "c" is for the insulated shield wires. 



Figure 33. The Maximum Magnetization Current for the Late-Time MHD-EMP 
Waveform when p = 100 fi-m. 
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Figure 33. The Maximum Magnetization Current for the Late Time MHD-EMP 
Waveform when p = 100 fl-m. 


63 






























b. Shield Wire Connected to Tower with Span of 0.40 km 

Figure 34. Time Delay to Transformer Saturation for the Late-Time MHD-EMP 
Waveform when p = 100 O-m. 
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c. Shield Wire Not Connected to Tower 

higvire 34. Time Delay to Transformer Saturation for the Late-Time MHD-EMP 
Waveform when and p = 100 Cl-m. 
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5.0 MHD-EMP EFFECTS ON THE POWER SYSTEM 


The MHD-EMP-induced currents on long lines can have an effect on the power 
system components connected to the lines [1]. The previous section has illustrated that 
these quasi-dc currents can cause transformer magnetization currents to increase. At the 
same time, the harmonic distortion of the 60-Hz power on the system increases, additional 
reactive power is required of the generator, and the power system relaxation time 
constant changes. Each of these phenomena can have an effect on the power system 
stability. 

5.1 Harmonic Generation 

Reference [13] describes a series of measurements made on a transmission-class 
system in which the harmonic content of the magnetization currents was measured with 
various levels of dc current injected into the transformer neutral. The transformers used 
in these tests were autotransformers. 

The first configuration tested consisted of a bank of single-phase 500/230-kV 
autotransformers having a 34.5 tertiary winding. These transformers were rated at 
120/230 MVA per phase, or 360/600 MVA for the complete three-phase unit. Figure 35 
presents the measured harmonic content on the primary winding of one of the 
transformers, which was unloaded on the secondary. These data are shown for various 
levels of dc injection into the transformer neutral. Because the transformer was unloaded, 
the fundamental component of the current at 60 Hz, is relatively small, and it is shown 
along with the higher harmonics. To be consistent with ref.[13], the component of the 
current at 60 Hz in this report is referred to as the first harmonic. The waveform 
component at 120 Hz is the second harmonic, and so forth. 

A second transformer which involved a smaller 230/115-kV autotransformer with a 
13.8 tertiary winding was also considered in ref. [13]. This transformer was rated at 
200/333 MVA, and had a three-legged core design. The phase conductor of the 
transformer which was measured had an electrical load attached. Figure 36 illustrates the 
harmonic content in the line current on the 230-kV side of this transformer for various 
levels of dc excitation. Because the resulting line current at the fundamental 60-Hz 
frequency was very high and off the scale, it is not shown on this plot. Only the second 
and and higher higher are presented. 
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Measurements of E 3 responses were also made on a mock-up of a three-phase 
480-V/ 12.47-kV/208-V power distribution system [14]. A 480-V/12.47-kV 300-kVA Y-Y 
step-up transformer was connected to a 12.47-kV/208-V 75-kVA Y-Y step-down 
transformer, and dc current was injected into the neutrals of the transformers. With the 
system energized by commercial 480-V power, measurements of the harmonic distortion 
were made. These results are presented in Figure 37, shown in the same format as the 
previous harmonic data. As in Figure 36, the 60-Hz component of the line power current 
was high (well over 2.0 A); consequently, it is not plotted in this figure. 

In these plots, it is clear that small dc currents on the order of 2-5 A can have a 
significant effect on the harmonic content of power distribution systems. Similarly, 
currents on the order of 50-100 A can affect the power transmission system. Figure 10 
indicates that for a 12-kV distribution system having a line length of only about 5 km, 
the normalized dc current I C /E Q induced by the MHD-EMP environment would be 
about 0.4 A-km/V. For an assumed peak MHD-EMP E-field strength of just 20 V/km, 
this results in a current of 8 A, which is above the levels of current observed to cause 
harmonic distortion in ref. [14]. Similarly, for a 500-kV power transmission line which is 
100 km long, Figure 9 indicates that I C /E Q a 40 A-km/V. Again for a peak MHD-EMP 
E-field strength of 20 V/km, this results in a current of 800 A, which is well above the 
current levels observed to produce harmonic distortion. Of course, a larger MHD-EMP 
field strength or a longer line will make these effects even more pronounced. Thus, it is 
evident that the MHD-EMP environment will cause serious harmonic distortion in both 
power transmission and distribution systems. 
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Figure 35. Measured Harmonic Content for Dc Injection in the Primary of an Unloaded 
500/230-kV Autotransformer, from Ref. [13]. 
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Figure 36. Measured Harmonic Content for Dc Injection in the Primary of a Loaded 
230/115-kV Autotransformer, from Ref. [13]. 
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Figure 37. Measured Harmonic Content for Dc Injection in the Neutral of a Three-Phase 
480-V/12.47-kV/208-V Power Distribution System, from Ref.[14]. 
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5.2 Reactive Power Demand Caused by MHD-EMP 


Reference [13] also reports the change in the reactive power flowing into the 
500/230-kV autotransformer under dc excitation conditions. Figure 38 shows this power 
in Mvar per phase, as a function of the dc level. Similarly, measurements of the reactive 
power on the distribution system in ref.[14] were also made, and these results are 
summarized in Figure 39. Note that the vertical scale in figure 39 is in Kvar/phase. 

Both of these increases in reactive power are significantly lower than the rated 
power levels of the transformers. Consequently, this increase in power will not be 
detrimental to the transformers. However, the power generation units will be required to 
produce this additional reactive power, and if the demand is too great, the generation 
facilities will not be able to provide sufficient power to maintain the proper operating 
voltage in the system. The voltage on the power system will consequently drop, and a 
power outage may ensue. 

5.3 Changes in System Time Constant 

The testing on the power distribution system in ref.[14] found that the 
characteristic relaxation time of the power system depends on the level of dc excitation. 
For a very small dc excitation, the time required for a transient perturbation in the 60-Hz 
current to decay (the time constant) can be on the order of several seconds. However, as 
the level of dc excitation increases, this time constant becomes less than one second. 
Figure 40 (from ref. [14]) illustrates the measured system time constant as a function of 
the dc injection. 

The results of Figure 40 indicate that the effect of the E 3 waveform on a power 
distribution system will occur within 1 second of the application of the field. The 
early-time part of the E 3 waveform shown in Figure 4a has typical rise and fall times that 
are longer than the typical system relaxation time of 0.5 second. The late-time E 3 in the 
figure is much slower than this relaxation time. Thus, for all practical purposes, the 
MHD-EMP response can be treated as a purely dc effect, neglecting the R/L time 
constants inherent in the transformers. 
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Figure 38. Reactive Power Demand for the 500/230-kV Transformer, from Ref. [13]. 



Dc Injection 

Figure 39. Reactive Power Demand for the 480-V/12.47-kV/208-V Power Distribution 
System, from Ref. [14]. 
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Figure 40. System Relaxation Time vs. DC Current Injection for the Three-Phase 
480-V/12.47-kV/208-V Power Distribution System, from Ref.[l4]. 
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The system relaxation time discussed above was for the distribution system with a 
nominal 60-Hz power load. Recent tests in Minnesota suggest that the electric load has a 
substantial effect on transformer saturation time constants [15]. This fact also has been 
verified with the calculational model described in this report. Specifically, Figures 4144 
illustrate the effects of the electric load on the saturation time constants. The model 
selected is a system with the following parameters: 


MHD-EMP Environment 
Soil Resistivity 
Line Length 
Span Length 

Tower Footing Resistance 



30 km 
0.16 km 
25 D 


The MHD-EMP waveform was assumed to be a step voltage to easily demonstrate the 
saturation time constant. Figure 41 illustrates the system response when the transformer 
is unloaded. Figures 42, 43, and 44 illustrate the system response when the transformer is 
loaded with 20%, 40%, and 100% load, respectively. The load connection is assumed to 
be wye-grounded. The saturation time constants are shown in Table 8. 


Table 8 


Saturation Time Constants for Loaded Transformer 


Waveform Type: 20 V/km Step Function 
Shield Connected to Towers 
Span Length: 0.16 km 
Ground Resistivity: 100 ft*m 


Transformer Load Saturation Time 


(in %) 


Constant (Sec) 


0 % 
20 % 
40% 
100 % 


48 

17 

11 

6 
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b. Transformer Magnetization Current 



c. Transformer Magnetic Flux 


Figure 41. System Response to MHD-EMP with Transformer Unloaded. 
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a. Line Zero-Sequence Current 




c. Transformer Magnetic Flux 

Figure 42. System Response to MHD-EMP with 20% Transformer Load. 
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b. Transformer Magnetization Current 



c. Transformer Magnetic Flux 

Figure 43. System Response to MHD-EMP with 40 % Transformer Load. 
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b. Transformer Magnetization Current 



c. Transformer Magnetic Flux 

Figure 44. System Response to MHD-EMP with 100% Transformer Load. 
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When the electric load was assumed to be delta-connected, no effect was observed 
on the saturation time constants. 

The observed results can be explained with a simplified model of the zero-sequence 
network. Specifically, one can think of the saturation time constant as determined by an 
equivalent R-L circuit model. In the equivalent model, the inductance is dominated by 
the transformer inductances, and the resistance is dominated by the transformer winding 
resistance and affected partly by the load resistance. The transformer magnetizing 
inductance can be thought of as being connected in parallel. For an autotransformer with 
tertiary winding, the influence of the electric load on the equivalent resistance of the 
zero-sequence model has been determined to be 

Zgq = ( z h * Z H ) + (2Z d + Z L ) (2Z t -t- Z L ) / (Z L + Z T + Z D ) (28) 

and for the wye-wye connected transformer with a tertiary winding, the resistance is 

Z eq = Z H + ( Z D + Z L^ Z T / ( Z L + Z T + Z D^ ( 29 ) 

In Eqs.(28) and (29) the following definitions are used: 

Zjj = Primary leakage impedance, 

Z^ = Secondary leakage impedance, 

Zrp = Tertiary leakage impedance, 

Zjy = Electric load impedance. 

The above formulae should be viewed only as approximations. They can provide 
the approximate increase of the equivalent resistance due to the electric load. As an 
example, assuming a 20% electric load for a transformer with leakage impedance of 0.001 
+ j 0.10 pu, the effect of the electric load is to double the equivalent resistance. Doubling 
the equivalent resistance will substantially decrease the saturation time constant due to 
the nonlinear characteristic of the magnetization inductance. 
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As in the case of the change in reactive power, these changes in the transformer 
time constants will not have a direct adverse effect on the power system components. 
System-wide, however, the overall effect of such changes in the operating parameters of 
the power system is not at all clear, and performing an assessment of these effects at this 
level is beyond the scope of this report. 

5.4 Observed Geomagnetic Storm Effects 

As mentioned earlier in this report, the MHD—EMP environment is similar to that 
of a geomagnetic storm. Consequently, it is useful to look to the several reported power 
system malfunctions which have been ascribed to these naturally-occurring storms. 
Doing this can provide an indication of possible MHD—EMP effects on the commercial 
power system. 

On October 28 and 29, 1991, a major geomagnetic storm occurred. This was given 
a K-Index of 9 from readings made in Boulder, Colorado, and in Loring, Maine. The 
storm began at 1540 universal time (1040 EST, 0940 CST, 0840 MST, 0740 PST) on 
October 28. A number of geomagnetic observatories across the US and Canada recorded 
the fluctuations of the geomagnetic field during this event and can provide data. Of 
particular interest are the data from the Canadian Geological survey of Canada [16], as 
the data are available at 10 second measuring intervals. Other data, such as those 
available from the U.S. Geological Survey, are typically available only with a one minute 
time resolution [17]. 

Figure 45 presents the measured North-South and East-West geomagnetic field at 
the Ottawa Magnetic Observatory, as provided by ref.[16]. These data are plotted as a 
continuous record for a three day period, starting on 10/27/91. The vertical axis is the 
absolute magnetic field (i.e., the static geomagnetic field plus a small time-varying 
component) in units of nano—teslas. The onset of the geomagnetic storm is defined to be 
at about 15:40 universal time on October 28, and this time is indicated on the plots. Note 
that the North—South B—field is the dominant component. 


The electric field on the earth's surface corresponding to these magnetic fields may 
be calculated by numerically evaluating a convolution integral, as described in Appendix 
A of ref. [ 1 ]. This assumes a simple, homogeneous conducting half-space as a model of the 
earth. Figure 46 presents the resulting calculated E-fields on the earth for an assumed 
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Figure 45. Measured Geomagnetic Field at Ottawa, Canada, from Ref. [16]. 
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electrical conductivity of <r — 0.0001 mhos/m. Part a of the figure shows the 
North-South E-field, which is derived from the East-West B-field. Part b illustrates 
the East-West E-field. It is interesting to note that although the North-South B-field 
is the larger of the magnetic field components, the East-West E-field is not the 
dominant electric field. There is a 14 V/km spike in the North-South E-field which 
caused several problems in power systems across the U.S. 

To observe in better detail the behavior of the earth-induced E-fields at the onset 
of the storm, Figure 47 presents the fields on an expanded time scale, in minutes after 
15:30 universal time on October 28. From this plot, it is clear that the 10 second time 
iesolution is not sufficiently small to adequately sample the fast—rising peaks in the 
E field components. These electric fields actually could have been larger in amplitude 
than what is indicated in this figure. 

This geomagnetic storm caused a number of power system problems which have 
been reported in [18]. These are reproduced below: 


• 2SSS N° rth ? ast Pow er Coordinating Council) - The Radisson-Sandy Pond 
Fm j w ne tnp , ped Monday morning (October 28) because filters on the line 
laded Very early Tuesday morning, the Madawaska HVDC tie between Quebec 
and New Brunswick tripped. Transformer problems were reported, but not 
confirmed as solar storm related. 

• MAAC (Mid—Atlantic Area Council) — PJM detected geomagnetic induced 
currents (GIG) due to the solar storms on several occasions. Late Monday night 
until early Tuesday morning (October 28-29), the GIC reached the Region's 
operating procedure "trigger level" and west—to—east electricity transfer limits 
were reduced. Public Service Electric & Gas Company also reduced the output 
ot its Salem nuclear generating plants. 

• MAIN (Mid—America Interconnected Network) — Transformers at the Point 
Beach nuclear units were "noisy" for a while, but the noise stopped when the 
system operators backed down on the var output of the generating units Hat 
this not worked, opening of a major 345 kV line from North appleton to Rocky 
Run would have been the next step to protect the transformers. MAIN also is 
checking to determine whether significant tile lie flow changes on Monday 
(October 28) and the tripping of two generating units were solar storm related. 

• MAPP (Mid-Continent Area Power Pool) - One utility in MAPP reported 
voltage fluctuations. 
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Figure 46. Computed Earth—Induced Electric Field for a = 0.0001 Mhos/m. 
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Figure 47. Computed E-Field Components for October 28, 1991. 


WSCC (Western Systems Coordinating Council) - BP A (Bonneville Power 
Administration) reported that on Monday (October 28) about 15 capacitor banks 
at seven substations (115 and 230 kV) tripped. The Blackwater HVDC tie 
between WSCC and SPP also tripped. Both British Columbia Hydro and 
that day* Ca lforma Ec * 1S0n reported voltage excursions lasting several minutes 


Furthermore, ref.[19] documents the following anomalies for the same geomagnetic storm: 


BPA (Bonneville Power Administration) - At 0737 PST October 28, 1991 shunt 
capacitors started tripping and continued to trip until 0739 PST at Chemwa 
Keeler, Tillamook, Bandon, McMinnville, Alvey, Cosmopolis, and Ling View 
substations. Operators reported "strange" transformer noises at Keeler and Peal 
Substations. Also, there was a report of transformer noise at Portland General 
Electric's Boardman Generation Plant. 


• SCE (Southern California Edison) - At 0737 PST October 28,1991, the 550 kV 
transmission system voltage dropped from 530—535 kV to 518—524 kV It 
recovered within three minutes. 


At 0801 PST October 28, 1991, the 550 kV transmission system voltage again 
dropped from 530—535 kV to 522—526 kV. It recovered within three minutes. 
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Also a fault recorder triggered on abnormal neutral current at the Serrano 
230/500 kV substation which is located southeast of Los Angeles. Off-line 
analysis indicated the neutral current magnitude was between 150 and 200 
amperes and principally third harmonic current with 32% 6th harmonic current. 
Phase A current evaluation revealed strong 2nd and 4th harmonic currents, 1.7% 
and 1.0% respectively, of the bank's rating. 

• NMPS (New Mexico Public Service) — At 0838 MST October 28, 1991, the 
ac—dc—ac back—to—back tie at Blackwater went off line. Preliminary indications 
are that it was tripped by an overvoltage relay set to trip at 115% voltage for one 
second. 

• NEPSCO (New England Public Service Co.) - At 1037 EST October 28, 1991, 
the New England Hydro Phase II dc tie with Hydro Quebec tripped. The 
problem caused by harmonics in the dc system appears to have occurred at the 
Radisson Terminal. NEPSCO was importing 930 MW at Sandy Pond at the 
time of the trip. The system was back in service within an hour. 

• PSE&G (Public Service Electric & Gas Company) — At 2257 EST October 28, 
1991, the Salem Unit was backed off to 80%. 

• APS (Allegheny Power System) - At 1050 EST October 28, 1991, Allegheny 

Power System SCADA armed the capacitor bank trip restrain and enabled the 

Meadowbrook transformer gas detector trip. No capacitor banks tripped in the 
minutes before the restrain was initiated. 

At 2259 EST October 28, 1991, the T4 transformers at Meadowbrook tripped on 
gas detection. 

At 2316 October 28, 1991, APS removed T2 at Meadowbrook per operating 
procedures. 

• WEPCO (Wisconsin Electric Power Co.) — Reported transformer growling at 
approximately 0930 CST October 28, 1991, at its Point Beach Power Plant. 
They backed off VAR generation during this period. 

• VEPCo (Virginia Electric and Power Co.) — At 1041 EST October 28 , 1991, 

230kV — 100MVAR capacitor bank at Chuckatuck tripped due to neutral 

unbalance. 

At 1041 EST October 28 , 1991, 230kV — 150MVAR capacitor bank at Dooms 
tripped due to neutral unbalance. 

At 1102 EST October 28 , 1991, 115kV — 25MVAR capacitor bank at Staunton 
tripped due to neutral unbalance. 
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From these observations of geomagnetic storm effects on the power system, it is 
reasonable to conclude that the power system would respond to a MHD-EMP event. 
Figure 48 presents an overlay of the earth—induced E—field for the geomagnetic storm on 
October 28, and a normalized composite MHD-EMP E-field, similar to that shown in 
Figure 4a. The MHD-EMP waveform has been normalized here to have the same peak 
amplitude of about 14 V/km. Note that there is a remarkable similarity in these two 
waveforms: an early—time spike, followed by a later-time component. 





6 7 8 9 10 11 12 

Time (Minutes) 


Figure 48. Comparison of MHD-EMP with Geomagnetic Storm E-field Waveforms. 
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6.0 SUMMARY AND CONCLUSIONS 


This report has discussed the interaction of electrical transmission and distribution 
systems with the MHD-EMP, or E 3 , environment. A simplified form of the earth-induced 
electric field has been postulated, and a number of different calculations of coupling to 
various line configurations have been performed. 

The time variations of the MHD-EMP electric field are sufficiently slow that an 
analysis of quasi-dc coupling to the lines is feasible. This analysis reduces the 
electromagnetic field coupling problem to that of a relatively simple electrical circuit. 
The power system reacts to the external MHD-EMP excitation with its own R/L time 
constant, which depends on the electrical configuration and parameters of the system 
being considered. For the MHD-EMP environment to be able to induce currents in the 
power system, a dc path to ground at both ends of a long length of electrical conductor is 
required. One way of obtaining this configuration is to use grounded-wye transformers at 
each end of the line. For an E 3 electric field strength of 20 V/km and the assumed line 
and transformer parameters, typical peak quasi-dc currents for a 500-km, 500-kV 
transmission line was found to be on the order of 1600 A. For a 50-km, 69-kV 
subtransmission line and for a 30-km, 12-kV distribution line, these MHD-EMP - induced 
currents are between 70 and 190 A, and between 40 and 110 A, respectively, depending on 
how well the transformers are grounded. Removing the ground connection of the 
grounded-wye transformer or using an ungrounded delta transformer in its place are two 
ways of eliminating this MHD-EMP current. 

For high-voltage transmission systems, we found in a sample calculation that times 
required to saturate loaded power transformers range from 6 to 11 seconds. For 
lower-voltage distribution-class systems, the saturation time can be on the order of a 
second. Response times for both of these systems are seen to decrease nonlinearly with an 
increase in the MHD-EMP field strength. 

The MHD-EMP environment has been shown to cause saturation in transformers 
in both transmission and distribution systems. Using measured transformer response data 
for a dc current injection into the neutrals, we have observed that quasi-dc current levels 
corresponding to E-field strengths as low as 20 V/km for typical line lengths can generate 
60-Hz harmonics. Thus, it is evident that the MHD-EMP environment will cause serious 
harmonic distortion in both power transmission and distribution systems. 
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We also note that there is an increase in the reactive power flow within the system 
when it is excited by the MHD-EMP environment. Specific data for the increase in 
reactive power demand for both transmission- and distribution-class transformers have 
been presented. This increase in reactive power arises due to the saturation of the 
transformers and can cause instabilities in power system operation. The severity and 
duration of these system instabilities is presently unknown. 

From this work, it is clear that MHD-EMP field strengths on the order of 20 
V/km can have measurable effects on transformer operation in the form of core 
saturation, harmonic generation, and an increase in reactive power demand. Geomagnetic 
storms with E—Field peaks on the order of just 10 to 15 V/km have been known to cause 
problems in the power system. Larger MHD-EMP E-field strengths will only increase the 
severity of these effects. The system-wide impact of MHD-EMP on a large power 
network in which several different transformers experience these effects simultaneously 
remains to be determined. 
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